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Abstract 

Localized surface plasmon resonances (LSPRs) have recently been identified in extremely 
diluted electron systems obtained by doping semiconductor quantum dots. Here we investigate 
the role that different surface effects, namely electronic spill-out and diffuse surface scattering, 
play in the optical properties of these ultra-low electron density nanosystems. Diffuse scatter¬ 
ing originates from imperfections or roughness at a microscopic scale on the surface. Using an 
electromagnetic theory that describes this mechanism in conjunction with a dielectric function 

*To whom correspondence should be addressed 

^Departamento de Fisica Teorica de la Materia Condensada C5 and Condensed Matter Physics Center (IFIMAC), 
Universidad Autonoma de Madrid, E-28049 Madrid, Spain 

^Centre of New Technologies, University of Warsaw, Banacha 2c, 02-097 Warsaw, Poland 
'ilDepartment of Applied Physics and Gothenburg Physics Centre, Chalmers University of Technology, SE-412 96 
Goteborg, Sweden 


1 




including the quantum size effect, we find that the LSPRs show an oscillatory behavior both in 
position and width for large particles and a strong blueshift in energy and an increased width 
for smaller radii, consistent with recent experimental results for photodoped ZnO nanocrystals. 
We thus show that the commonly ignored process of diffuse surface scattering is a more impor¬ 
tant mechanism affecting the plasmonic properties of ultra-low electron density nanoparticles 
than the spill-out effect. 
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Surface plasmons in metals have been studied for a long time.^i^ However, in the last few 
years, the behavior of plasmon resonanees in nanosystems has beeome a hot topie due to numer¬ 
ous applieations,® having been made possible by eonsiderable progress in both ealeulating and 
measuring the optieal properties of nanopartieles. A partieularly interesting example eoneems the 
energy shift experieneed by the loealized surfaee plasmon resonanees (LSPRs) of nanospheres as 
a funetion of the size, as well as the eorresponding modifieation in the plasmon width.^IH! 

Reeently, it was diseovered that doped semieonduetor quantum dots with a high eoneentration 
of earriers in the eonduetion band show, when illuminated, absorption resonanees attributed to the 
exeitation of surfaee plasmons Not only are doped semieonduetor nanoerystals emerging as 
a new elass of plasmonie struetures for their remarkable properties of tuning and switehing on/off 
the plasmon energyl^ The density of eleetrons in sueh nanoerystals is low, typieally two to 
three orders of magnitude lower than in typieal metals. In terms of the elassie density parameter 
A = ( 4 ^)^’ being the earner density, whieh in units of the Bohr radius is 2-6 for eommon 
metals, in doped semieonduetors is on the order of 25. Sueh a diluted eleetron gas has its 
absorption band in the infrared rather than in the visible part of the eleetromagnetie speetrum. 

In previous publieations^^^^ we analyzed how the so-ealled surfaee spill-out effeet modifies 
the position and width of plasmonie resonanees in Ag and Au nanospheres of radii smaller than 
10 nm embedded in different dieleetrie matriees. This effeet originates from the faet that the elee- 
tronie density of typieal high eleetron density metals does not have a sharp diseontinuity at the 
geometrieal surfaee, as assumed by the elassieal theory. Rather, it is a eontinuous funetion beeause 
the surfaee potential barrier eonfining the eleetrons is finite and smooth. These studies omitted 
another surfaee meehanism whieh affeets loeal surfaee plasmon resonanees in any geometry, and 
is assoeiated with the existenee of imperfeetions or roughness at a mieroseopie seale on the sur¬ 
faee. As a eonsequenee, some of the metal eleetrons arriving at the surfaee are refleeted baek in 
any direetion.^^^J^^ The origin of the phenomenon ean be better understood in the ease of planar 
surfaees. For an ideally flat surfaee, the translational invarianee of spaee in the direetion parallel to 
the surfaee implies eonservation of the parallel momentum and, therefore, eonservation of the par- 
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allel component of the eleetrie eurrent (speeular seattering). However, this restrietion is lifted for 
non-ideal surfaees, where it is possible to have diffuse seattering, and, eonsequently, fluetuations 
in the parallel eurrent. This meehanism affeets any physieal phenomenon taking plaee at a surfaee 
and, in partieular, will lead to energy broadening of the LSPRs sinee it breaks eoherenee of single 
seattering events. For planar surfaees or spheres with typieal metallie eleetronie densities, the spill- 
out meehanism dominates over the diffuse surfaee seattering one at the frequeneies of the surfaee 
plasmon modes, this being the reason why the seeond meehanism is generally negleeted.^^^ How¬ 
ever, imperfeetions in the shape and morphology of ultra small partieles are frequent and eonstitute 
an important souree for diffuse surfaee seattering, partieularly for the low-density eleetron systems 
that we will study in this work. Similar behavior ean be produeed by adsorbed moleeules.^^^l^ll 

The purpose of this artiele is to investigate the role that spill-out and diffuse surfaee seattering 
play in the optieal properties of ultra-low eleetron density nanosystems. Using an eleetromag- 
netie theory that deseribes diffuse surfaee seattering in eonjunetion with a dieleetrie funetion that 
ineludes the effeet of quantization due to the small size of the partieles, we find that the LSPRs 
show an oseillatory behavior both in position and width for partieles larger than ea. 3 nm in radius 
and a strong blueshift in energy and an inereased width for smaller radii, eonsistent with reeent 
experimental results.^ We thus show that the eommonly ignored proeess of diffuse surfaee seat¬ 
tering is an important meehanism affeeting the plasmonie properties of these systems. Based on 
our theory, we ean also extraet an effeetive length eharaeterizing diffuse surfaee seattering whieh 
ean be ineluded in an effeetive dieleetrie funetion and used in simpler models. 

Surfaee seattering meehanisms affeeting the eleetromagnetie response of metal surfaees ean be 
deseribed by means of two effeetive eomplex lengths, d^{G)) and (m), where (O is the frequeney. 
These are surfaee response funetions assoeiated with ehanges aeross the surfaee of the normal 
eomponent of the eleetrie field veetor and of the parallel eomponent of the displaeement veetor 
with respeet to their elassieal eounterparts, respeetively.^^The length d^{(0) ean be related to the 
eleetronie eharge density, 5p, indueed at the surfaee by any external perturbation and deseribes 
the spill-out effeet, while diffuse surfaee seattering, whieh provokes fluetuations in the eurrent 
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parallel to the metal surfaee, is deseribed by (<ij|((0) = 0 for speeular scattering).!^ For the 

case of a sphere of radius R, the corresponding lengths are denoted as dr{G),R) and dQ{G),R) and, 
by considering the mentioned changes in Er and Dq across the sphere surface respectively, it was 
shown!^^!^ that, in the quasi-static limit -C 1, the polarizability of the sphere can be written as 


a{(o) =R^ 




£{CO) +2£m + 2{£{(0) - £n 


dAa,R) I ' 


( 1 ) 


R 


Then, in the same quasi-static limit, the optical absorption cross section can be calculated as 


o{(o)=AK^^/^\m[a{(o)], (2) 

where c is the speed of light. In eqs. (??) and (??), e((0) is the classical, local permittivity of 
the metal and is a frequency independent permittivity of the surrounding medium. Note that if 
surface effects are neglected {dr = 0 and = 0 in eq. (??)), the classical results of the Mie theory 
are recovered. 

The role that these two lengths play in the effective relaxation time of the surface plasmons 
in small spheres with typical electron densities characterized by —2, 3 and 4, was analyzed 
in Ref.!^ At these densities the spill-out contribution to the surface plasmon relaxation time is 
much larger than the diffuse surface scattering contribution, same as for planar surfaces. However, 
as the electron density decreases, the spill-out contribution decreases quickly while the diffuse 
surface scattering one remains nearly constant (we refer the reader to Fig. 2 of Ref.Therefore, 
one expects the relative importance of both contributions to be reversed for low enough electron 
densities. 

These contributions can be estimated rather easily from eq. (??) and the results of Ref.!^(see 
Supporting Information, SI, for a detailed derivation). Assuming £„, = 1 and £{co) = 1 — (^)^, (Op 
being the bulk plasma frequency, we have the classical Mie LSPR frequency of the sphere cof = 
(Op{\/3. From eqs. (??) and (??) one can obtain its full-width at half-maximum (see eq. (S6)). 
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Substituting dr and dg by their counterparts for a planar surface {d^{(0) and d\\ (ft)), respectively), 
a good approximation for large R, we find that for a perfectly diffuse surface reflecting electrons at 
random, the ratio of spill-out to diffuse surface scattering contributions can be estimated as 



(3) 


which is cast in terms of the effective electron mass m* (in units of the electron mass nig) and the 
one-electron radius (in units of the Bohr radius qq), with qq = 


Equation (??) shows that the spill-out contribution decreases for effective masses smaller than 
the electron mass, which is the case for the conduction electrons in doped semiconductors. Also, 
the ratio of the contributions of the two mechanisms to the plasmonic width decreases with in¬ 
creasing Vs as lm[—d^{ 0 )f)]/y/Fs. Taking the values of Im[—for = 2 - 5 reported in 


Refand m* = 1, one obtains values of decreasing from 2.5 for = 2, to 0.5 for rs — 5. 


This estimate, as well as a simple analysis in the SI, show that the effects of non-specular scatter¬ 
ing become more important as the electron density is lowered and will largely overcome these of 
electronic spill-out for ~ 25. 

Having demonstrated the importance of accounting for diffuse surface scattering, we now pro¬ 
ceed to analyze these effects in a detailed manner. A microscopic description of this effect is, 
however, beyond the scope of the present work. Instead we use a phenomenological theory first 
designed for planar surfaces^^^^^ and then extended to spheres.!^ The theory embeds the real 
sphere in an infinite, fictitious medium having exactly the same dielectric functions. Then, a con¬ 
stitutive relation giving the polarization Pf (r, to) due to the free charges inside the real sphere, is 
written as 



(4) 


where is the electric field vector in the infinite medium, which, for |r| < R is the actual electric 
field inside the sphere, Cq is the permittivity of free space, £ is the dielectric tensor of the medium. 
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Cool accounts for interband transitions of bound electrons, and the integral extends to the whole 
space. In this integral, the region of space |r'| < R describes the excitations produced at a point 
r' inside the sphere that propagate directly to the point r while the fictitious region |r'| > R sim¬ 
ulates these excitations that arrive at r after being reflected at the surface. Therefore, the surface 
properties are mimicked by the values of the electric field in the fictitious region of the infinite 
medium. In this work, we want the sphere surface to reflect electrons completely at random and 
this means that, on the average, no excitation will arrive to r coming from the surface. In eq. (??), 
this is thus equivalent to making the electric field to be zero outside the sphere. Hence our 
problem consists of constructing an electric field of the form 

E^D/!ej'e(r, w) for |r|<7? 

I 0 for |r| > R 

satisfying the Maxwell equations. In general, this is not possible without adding fictitious charges/currents 
to the ficticious region of the medium. Here we will use a different, equivalent, procedure. 

From the point of view of an electromagnetic theoretical formulation of both kinds of surface 
effects discussed above, one should note that the radial component deviates from its classical 
counterpart because of the existence of longitudinal modes localized in the surface region, which 
are not present in the classical formulation. Consequently, the spill-out effect can only be analyzed 
if a non-local longitudinal dielectric function, that is, dependent on the spatial coordinates as well 
as on the frequency, is included in the theory. In a similar way, a description of the diffuse surface 
scattering mechanism requires the use of a non-local transverse dielectric function, which, as we 
will see below, allows to include in the theory the excitation of transverse electron-holes pairs in 
addition to the classical polariton mode propagating with wave vector k; = “ a/£(w). Since our 
previous analysis shows that we can expect the spill-out effects to be small, we will neglect them 
completely by making the longitudinal dielectric function of the sphere a local one. However, the 
transverse dielectric function has to be non-local. In our approximation £ only depends on spatial 
coordinates through the difference r — r', and it is convenient to Fourier-transform the permittiv- 
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ity to momentum space. We will use here the simplest possible form for a non-local transverse 
dielectric function 


er(k,(o) = 


co^-A'^ + icoYh - ’ 


(5) 


where k is a wave vector, |k| = and Pt is a constant proportional to the Fermi velocity Vf. This 
form of £t is obtained when the whole spectrum of electron-hole pairs is substituted by just a 
single pair. Thus the pole of eq. (??) describes a transverse electron-hole pair dispersing in energy 
linearly with the wave vector. Quantum size effects (QSE) are included in the theory by means of 
the energy gap A = cOp^ti^with Rq = = (3;r^ne)3 being the Fermi wave vector. In 

this model, the strength of the QSE with particle size is proportional to Rq/R. 

The longitudinal and transverse dielectric functions have to be equal in the k = 0 limit and 
therefore a good approximation for £p is 


£l{G}) = 


(On 


( 6 ) 


It is convenient to solve for in a basis set of spherical vectors functions: (A:, r), mp/„j (fc, r) 

and \\.pim{k^r) represented by their quantum numbers, Z,m and the parity p (even or odd) with re¬ 
spect to the azimuthal angle. Here I and m only take positive integer values. The electric field 
vector of a plane wave of wave vector f incident onto the sphere along the x direction, 
can be expressed with components having simplyfying the problem. Then, the field in 

the extended medium can be written in general as 


(7) 

The field inside the metal sphere is a linear combination of the different normal modes that can 
be excited according to the proposed dielectric functions. In our case, there are no longitudinal 
modes because £i is a local dielectric function and there are two transverse modes (index j = 1,2): 
the polariton-like mode of wave vector T 2 and one transverse electron-hole pair of wave vector T\ 
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(see SI for details). Thus, the field inside the extended medium ean be written as 



for r<R 

0 for r > i? 


(8) 


in) 

At this point, the eleetric field inside the sphere eontains four eonstants, namely, E) / and 
0 = 15 2) and the electric field outside the sphere contains two constants for the scattered 
amplitudes of the even and odd components. However, we have four matching conditions for the 
electromagnetic field at the sphere surface. The missing two equations can be obtained by imposing 
that the free-polarization Vf calculated from the definition of eq. (??) with the field Ej'^(r) fulfills 
the wave equation V x (V x E) — (® )^ecx,E = inside the sphere (details in the SI). 

The electric field in the medium outside the sphere, has an odd component given by 


°° 2 / +1 

Eo(r) = J^/'^^y^[mo/i(fc„,r)+a/m^;i(fc„,r)], 
and the even component is 

°° 2 / +1 

Solving for a/ and bi allows us to calculate the absorption cross section 


(9) 


( 10 ) 


~J^(2/+l)Re[a/ + h/], (11) 

l=\ 

as scattering can be neglected in all the results we present due to the small size of our particles. In 
the following calculations we use the values of the parameters appropriate for ZnO nanocrystals in 
toluene: = 3.72, £„ = 2.5, m* = 0.28 and jh =0.1 eV.^With these values, we use the lower 

limit of the experimental electron density, Ug — f x 10^° cm^^, to get the experimental surface 
plasmon energies at the largest radii. Then, = 25.3, the Fermi velocity vp = 0.59 x 10^ m s“^ 
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“<J = 0-70 eV. 

To show clearly the effects a diffuse scattering surface produces in the plasmonic resonances, 
we first perform a calculation of the absorption cross section neglecting quantum size effects (A = 
0) and with Pr = vp/y/S. This value of Pt is given by the long wave length limit of the Lindhard 
transverse dielectric function. The results are depicted in Fig. [^for spheres of radii R = 4, 3 and 
2 nm, compared to the classical results of the Mie theory.^ We can appreciate the widening of 
the resonances and the appearance of small oscillations on their right side which are caused by the 
additional propagating transverse mode. The resonance strongly blue-shifts in energy for = 2 
nm. 



Figure 1: Effects of diffuse surface scattering as seen in absorption. Absorption cross-section 
normalized to surface area as a function of energy for R = 4, 3 and 2 nm from top to bottom. 
Quantum size effects are neglected (A = 0). Dashed lines: classical Mie theory; continuous lines: 
present theory for a totally diffuse surface. Notice the asymmetric widening of the resonance and 
the strong blue-shift for R = 2 nm. The inset shows a magnification of the plot in the energy region 
above the peak (semi-log scale) to highlight the oscillations. 


From our calculated cross section we extract the position, cOr, and the width of the resonance, 
r^, defined as the full width at half maximum. Pt and A are taken as adjustable parameters given 
that they have to be increased from the values Pt = vp/y/S and A = cOp^, Rq = = 0.56 

nm, by less than 30% in order to reproduce the experimental results of Ref.^ A reason for this is 
that the value of Pt is justified for a large system. Moreover, other values of Pt can be found in 
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the literature sinee the form of eq. (??) substitutes the whole speetrum of transverse eleetron-hole 
pairs by a single pair. With respeet to A, a loeal dieleetrie funetion for few eleetrons in a box has 
been ealeulated.l^ It has the form of eq. (??) with At,ox = where Ng is the number 

of eleetrons and a is the side of the box. Our frequeney gap A seales with Ne,m* and partiele size 
exaetly with the same powers and a prefaetor that differs by 15% if the eleetron density is the same 
in both models. Sinee the prefaetor has to depend on the shape of the nanopartiele, a ehange of ea. 
20% is justified. 



Figure 2: The energy of the plasmonie resonanee versus the sphere radius R for four groups of the 
values of the parameters jSy and Rq. Dashed line: Pt = vf/VS and Ro = 0 (no QSE); dot-dashed 
line Pt = vp/y/S andi?o = 0-56 nm; eontinuous line: fir = l-'ivp/y/S and7?o = 0.56 nm and dotted 
line: Pt = 1.2vf / y/5 and = 1 -2 x 0.56 nm. The experimental values of Ref.^ are represented 
by the blaek squares. By ineluding QSE and slightly inereasing the parameters involved from the 
basie theory there is a good aeeount for the experimental findings. 


Eigure shows cor as a funetion of the sphere radius. The dashed line displays the surfaee 
plasmon energy for A = 0 and Pt = Vf I y/5; it presents small oseillations with R for radii larger 
that ~ 3 nm, and then quiekly blue shifts for smaller radii. When we inelude the quantum size effeet 
using Rq = 0.56 nm (dot-dashed line) we observe the additional blue-shift eaused by this effeet, 
whieh is more pronouneed for the smaller radii, but the oseillations remain. If we now inerease Pt 
by 30% keeping the same value of Rq (solid line) we further inerease the blue-shift at small radii 
and also ehange the oseillations at large radii. When we inerease both, Pt and Rq by 20% we obtain 
the dotted line whieh does not differ mueh from the previous ease. In both eases the experimental 
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results^ can be reproduced remarkably well. Thus a diffuse surface produces a non-negligible 
blue-shift of the plasmonic resonances in small spheres with a low density of electrons, a result 
which is not obtained in the simple theory of Refwhere they would be essentially unshifted. 



Figure 3: The width of the plasmonic resonace versus \/R for the the same groups of the values of 
the parameters as in Fig. dashed line: Pt = vf/ V5 and Rq = 0; dot-dashed line Pt = vf j a/5 and 
Ro = 0-56 nm; continuous line: Pt = l-'ivp/ a/5 and Ro = 0.56 nm and dotted line: Pj = I ■2vf/ a/5 
and Rq = \ .2x 0.56 nm. The experimental values of Ref.^^ are represented by the black squares. 
While in the simple model ofl^the plasmon width would increase linearly with 1/R for all radii, 
the detailed model for diffuse surface scattering produces this behavior only at large radii. The 
width actually oscillates at intermediate radii and then increases fast at short radii. Notice that the 
plasmon width is much more sensitive to the values chosen for Pt than the plasmon energy. The 
experimental width can be reproduced by the same set of parameters that reproduce the plasmon 
resonance. 


Figurej^depicts the width of the resonance as a function of 1/R, for the same set of values of Pt 
and Rq as in Fig. Note that, while in the pure classical theory the width of the resonance would 
be equal to the bulk value = 0.1 eV (bottom line of Fig. [^), the surface scattering mechanism 
increases the plasmon damping linearly in 1/R for large radii. Then, for 6 nm > R > 3 nm the 
width oscillates between 0.13 eV and 0.16 eV and finally, for radii smaller than ca. 3-2.5 nm the 
plasmon width increases again almost linearly with l/R with a higher slope which also increases 
with Pt- The experimental results of Ref.^ are shown by the full squares. The width is better 
reproduced by our calculation with the largest value of Pt- Note also that the plasmon width is 
much more sensitive to the values of Pt than the plasmon energy. 
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We now obtain the length de{co) from our eleetromagnetie theory. This ean be done analytieally 
given the simplieity of our transverse dieleetrie funetion and that we only need to eonsider the Z = 1 
eomponent of the fields having even symmetry. Using the small-argument asymptotie values of the 
spherieal Bessel funetions appearing in the eleetrie fields of the even mode (see eq. (S36)) we get 
de as 

= -i(e-e„)(iR]h\'\iR)j,(tR), (12) 

K 

where t is the pole of the transverse dieleetrie permittivity, eq. (??), and e = e 7 ’(fc = 0, (o). The 
origin of the oseillations in (Or and Tr shown in Figs, [^andj^is elearly related to the Bessel 
funetions appearing in eq. (??). Sinee tRc:::^ ^ (eq. (S27)), then \tR\ 3> 1 for the values of R used 
here and large-argument asymptotie expansions of the Bessel funetions ean be performed yielding 

de{<o,R)^i - ^^ll+e'^% (13) 

Comparing this equation with the eorresponding d\\ for the planar surfaee, eq. (S13), whieh is 
used in simple models, we see that the eorreetion to ^Z[| for a real sphere is the faetor 1 so 

a large Im[t] is required to get rid of the oseillatory faetor to have the same limit, assuming A is 
unimportant. Then the erossover from large values of R, where the sphere ean still be deseribed 
by the planar surfaee, to small values of R, is given by 2R ~ (Im[t])^^ The physieal meaning 
of this relation is elear. The deeay length of the transverse eleetron-hole pair mode is (Im[t])^^ 
When this length beeomes shorter than the sphere diameter 2R, the mode only feels one surfaee 
and eonsequently, the limit of the planar surfaee is reaehed. 

Figurej^shows the real and the imaginary parts of de as a funetion of R eompared to their large- 
argument asymptotie values given by eq. (??), for a typieal value of (Ot = 0.3 eV, jir = vp/y/S 
and Rq = 0 (no QSE). The non-zero values of Re[Je] at large values of R are the eonsequenee 
of having ^ = y in this ealeulation, whieh, however, produees negligible shifts of the surfaee 
plasmon energies shown in Fig. Note that, even if QSE are negleeted, de depends strongly on R 
for small radii. One ean note that the expansion of eq. (??) is a good approximation to de for radii 
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Figure 4: (a) The real part and (b) the imaginary part of de versus R for Pt =vf/ VS, Ro = 0 and 
COr = 0.3 eV. Red solid line: results of the full theory, eq. (??), and dashed line: results of the 
asymtotic expansion eq. (??). The continuous horizontal lines are the values of Re[J[|((0f)] and 
Im[J||((p,)], being good asymptotes for de. 

larger than ca. 2 nm. 

Surface scattering is usually included in simple theories by means of an effective damping rate 
Y= Ye +A^, A being a constant of the order of 1. This procedure produces a linear increase of the 
LSPR width with 1 /R and, for low electron density systems, a continuous red shift of its spectral 
position with decreasing While this may be a reasonable approximation for large radii, when 
Re[de]/R 1 and the main effect of surface scattering is to increase the plasmon width linearly 
with l/R while the plasmon energy presents small oscillations around its classical value (Figs. 
and [^respectively), an effective damping rate theory cannot capture the nature of diffuse surface 
scattering at short radii. As already noted, is a complex quantity whose real part cannot 

in general be neglected. A better approximation can be obtained by defining an effective size- 
dependent dielectric function for a sphere containing diffuse surface scattering, e(( 0 ,R), via the 
Rayleigh polarizability (see eq. (S45)) as 

£(m,R) = (14) 

and de{(0,R) expressed by eq. (??). Such a definition can be useful in the analysis of problems 


14 









involving the dielectric function of the nanosphere. 

In this article, we have investigated how different surface effects, namely electronic spill-out 
and diffuse scattering, impact in the optical properties of ultra-low electron density nanospheres 
which are experimentally accessible by electron doping of semiconductor quantum dots. We first 
estimated the relative contributions of spill-out to diffuse surface scattering, showing that the ratio 
decreases with effective electron mass and with increasing r^. Next, we used a more elaborated 
theoretical model for including the diffuse surface scattering mechanism into an electromagnetic 
theory. Using this model, we calculated the extinction cross-section and found that the plasmonic 
resonances show an oscillatory behavior both in position and width for sizes larger than ca. 3 
nm in radius and a strong blueshift in energy and an increased width for smaller radii. When 
our model is used in conjunction with a dielectric function that includes effects of quantization 
due to the small size of the particles, we are able to reproduce the position and widths of the 
measured resonances!^ using reasonable values of the two basic parameters, Pt and A, involved 
in the theory. From our formalism, we can extract the length de{co,R) that allows to include the 
effects of a diffuse surface into an effective size-dependent dielectric function, which can be useful 
for future work. We conclude that diffuse scattering at the surface of ultra-low electron density 
nanocrystals is an important mechanism affecting their plasmonic properties and has to be taken 
into account when designing plasmonic devices based on doped semiconductor nano structures. 
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